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ABSTRACT 

This paper describes the capabilities and operational utility of a version of the Mesoscale Atmospheric Simula- 
tion System (MASS) that has been developed to support operational weather forecasting at the Kennedy Space Cen- 
ter (KSC) and Cape Canaveral Air Station (CCAS). The implementation of local, mesoscale modeling systems at 
KSC/CCAS is designed to provide detailed short-range (< 24 h) forecasts of winds, clouds, and hazardous weather 
such as thunderstorms. Short-range forecasting is a challenge for daily operations, and manned and unmanned launches 
since KSC/CCAS is located in central Florida where the weather during the warm season is dominated by mesoscale 
circulations like the sea breeze. 

For this application, MASS has been modified to run on a Stardent 3000 workstation. Workstation-based, real- 
time numerical modeling requires a compromise between the requirement to run the system fast enough so that the 
output can be used before expiration balanced against the desire to improve the simulations by increasing resolution 
and using more detailed physical parameterizations. It is now feasible to run high-resolution mesoscale models such 
as MASS on local workstations to provide timely forecasts at a fraction of the cost required to run these models on 
mainframe supercomputers. 

MASS has been running in the Applied Meteorology Unit (AMU) at KSC/CCAS since January 1994 for the pur- 
pose of system evaluation. In March 1995, the AMU began sending real-time MASS output to the forecasters and 
meteorologists at CCAS, Spaceflight Meteorology Group (Johnson Space Center, Houston, Texas), and the National 
Weather Service (Melbourne, Florida). However, MASS is not yet an operational system. The final decision whether 
to transition MASS for operational use will depend on a combination of forecaster feedback, the AMU’s final evalu- 
ation results, and the life-cycle costs of the operational system. 



1. Introduction 

The National Aeronautics and Space Administra- 
tion (NASA) and United States Air Force (USAF) 
have been conducting ground and spaceflight opera- 
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tions at the Kennedy Space Center (KSC) and East- 
ern Range at Cape Canaveral Air Station (CCAS) 
since the early 1960s. Weather support to operations 
at KSC/CCAS requires detailed forecasts of winds, 
clouds, ceilings, fog, and hazardous weather such 
as thunderstorms. Forecasting these parameters for 
KSC/CCAS is a challenging task since the central 
Florida facilities are located in an environment where 
there is an absence of significant large-scale dynami- 
cal forcing during much of the year. Under these con- 
ditions, regional and local factors such as land-water 
boundaries, land use, vegetation type/density, and soil 
moisture play a dominant role in determining the 
short-term evolution of weather conditions (Pielke 
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et al. 1991; Xian and Pielke 1991; McCumber and 
Pielke 1981). Hence, guidance from current genera- 
tion global and regional models is of limited value for 
these forecasting problems. 

The implementation of mesoscale modeling sys- 
tems locally at KSC/CCAS is ultimately intended to 
provide accurate forecasts of specific thunderstorm- 
related phenomena such as lightning, precipitation, 
and high winds. These forecasts are important for re- 
ducing downtime due to false weather advisories and 
alerts and minimizing the impact on personnel and 
equipment due to hazardous weather events occurring 
without warning. Improved forecast reliability may 
also permit safe relaxation of weather-related launch 
commit criteria for manned and unmanned space 
launches and flight rules for shuttle landings. 

To meet the forecasting needs at KSC/CCAS, 
NASA funded Mesoscale Environmental Simulations 
and Operations (MESO), Inc. to develop a version of 
the Mesoscale Atmospheric Simulation System 
(MASS) configured specifically for short-range fore- 
casting in the vicinity of KSC/CCAS. In this configu- 
ration, the model is run at a horizontal resolution of 
1 1 km with physical parameterizations for precipita- 
tion, radiation, and surface hydrology physics that are 
designed to predict convection induced by local varia- 
tions in surface heat and moisture fluxes, and cloud 
shading. 

MASS is one of several mesoscale modeling sys- 
tems that has been used for both basic and applied 
research including model development (e.g., Kaplan 
et al. 1982; Uccellini et al. 1987) and detailed case 
studies (e.g., Kaplan et al. 1985; Zack and Kaplan 
1987; Manobianco et al. 1991; Manobianco et al. 
1994) for more than 10 years. Other mesoscale mod- 
els that have been used for numerous phenomenologi- 
cal and modeling studies include the Pennsylvania 
State University/National Center for Atmospheric 
Research (PSU/NCAR) mesoscale model (e.g., Kuo 
et al. 1992; Doyle and Warner 1993a), the Regional 
Atmospheric Modeling System (RAMS) (e.g., Pielke 
1974; Cram et al. 1992), and the Naval Research 
Laboratory’s limited area model (e.g.. Holt et al. 
1990; Chang et al. 1993). 

State-of-the-art mesoscale modeling systems typi- 
cally contain detailed physical parameterizations and 
are run at very high horizontal and vertical resolu- 
tions. As a result, the models require large memory 
and processing capabilities, and until recently, could 
only be run on the fastest supercomputing platforms 
such as the CRAY-YMP. However, the development 


of computer workstations during the past five years 
with sufficient memory and processing speed has 
permitted mesoscale models to generate real-time 
forecasts at a fraction of the financial cost that would 
be required to run these models on mainframe 
supercomputers (Buzbee 1993). 

The version of the MASS model developed to sup- 
port operational weather forecasting at KSC/CCAS 
is designed to run in real time on high performance 
workstations. There are other applications where 
mesoscale model forecasts are performed in real time 
on mainframe or workstation-class computers. For 
example, the PSU/NCAR mesoscale model is run 
routinely on an IBM 3090 to support research, teach- 
ing, and public service in the Department of Meteo- 
rology at PSU (Warner and Seaman 1990). The 
RAMS model has also been run in real time on a 
Stardent 3040 workstation at the Colorado State 
University (CSU) to generate mesoscale forecasts 
during the 1991-1992 winter season over the Colo- 
rado Rocky Mountains (Cotton et al. 1994). In addi- 
tion, the Forecast Systems Laboratory (FSL) runs 
RAMS on an IBM RISC 6000/Model 580 worksta- 
tion to produce real-time forecasts that are displayed 
at FSL’s in-house daily weather briefings (Snook 
et al. 1995). Table 1 highlights the basic attributes of 
the real-time mesoscale modeling systems at PSU, 
CSU, and FSL. 

At KSC/CCAS, RAMS is run without precip- 
itation or cloud physics in real time on an IBM 
RISC 6000/Model 550 workstation as part of the 
Emergency Response Dose Assessment System 
(ERDAS). In ERDAS, the resulting wind fields from 
RAMS are used to drive dispersion models that could 
provide emergency response guidance to operations 
at KSC/CCAS in case of a hazardous material release 
or an aborted vehicle launch (Tremback et al. 1994). 
ERDAS is currently being evaluated for operations 
at KSC/CCAS. Additionally, a parallelized version of 
RAMS with precipitation and cloud physics is being 
developed by the ASTER Division of Mission Re- 
search Corporation to run on a cluster of workstations. 
The Parallelized RAMS Operational Weather Simu- 
lation System (PROWESS) is designed specifically 
for localized and sea-breeze thunderstorm forecast- 
ing at KSC/CCAS (Lyons et al. 1994). 

The MASS and RAMS models are currently run 
on workstations located in the Applied Meteorology 
Unit (AMU) at KSC/CCAS. The AMU was formed 
in September 1991 by a triagency Memorandum of 
Understanding among NASA, the USAF, and the Na- 
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Table 1. Basic attributes of selected real-time mesoscale modeling systems. See text for discussion and references listed for 
additional information. 

Model 

(reference) 

PSU/NCAR 

(Warner and Seaman 1990) 

RAMS/CSU 
(Cotton et al. 1994) 

RAMS/FSL 
(Snook et al. 1995) 

Horizontal grid 
dimensions 
(grid resolution) 

41 x41 (90 km) 
43 x 43 (30 km) 

36x28 <100 km) 
34 x 30 <25 km) 

61 x 61 (10 km) 

Number of 
vertical levels 

15 

25 

25 

Boundary-layer 

physics 

Blackadar 

(Zhang and Anthes 1982) 

Smagorinsky deformation 
K closure 

(Smagorinsky 1963) 

Smagorinsky deformation 
K closure 

(Smagorinsky 1963) 

Grid-scale 

precipitation 

physics 

Condense water vapor in 
excess of supersaturation 
to form precipitation 

Condense water vapor in 
excess of supersaturation; 
dump a portion as precipitation 
(Flatau et al. 1989) 

Condense water vapor in 
excess of supersaturation 
to cloud water 

No precipitation 

Subgrid-scale 
precipitation physics 
(cu mulus 
parameterization) 

Kuo type 

(Kuo 1974; Anthes 1977) 

None 

None 

Radiation physics 

Shortwave and longwave 
radiative surface fluxes adjusted 
based on clouds specified 
as a function of predicted 
relative humidity 

Radiative effects of water vapor, 
ozone, and carbon dioxide 
(Mahrer and Pielke 1977) 

Radiative effects of water 
vapor and clouds 
(Mahrer and Pielke 1977; 
Thompson 1993) 


tional Weather Service (NWS). The AMU’s mission 
is to evaluate and transfer new technology into the 
real-time weather support environment for the USAF 
Range Weather Operations (RWO) at CCAS and the 
National Oceanic and Atmospheric Administration's 
(NOAA’s) Spaceflight Meteorology Group (SMG) at 
the Johnson Space Center. The RWO and SMG are 
the two agencies that provide direct weather support 
to NASA. As part of a coordinated technology tran- 
sition plan for mesoscale numerical modeling, the 
AMU is currently evaluating MASS and ERDAS, and 
will be evaluating PROWESS. Based on the results of 
the evaluations, the AMU will recommend and/or de- 
velop modifications to these systems to make them suit- 
able for operational use and work with the Eastern 
Range to generate documentation and training nec- 
essary for certified Range systems. 

The objective ot this paper is to describe the gen- 
eral capabilities and operational utility of the version 
of MASS, which has been designed specifically for 
short-range forecasting at KSC/CCAS. Section 2 pro- 


vides an overview of the system and discusses the sta- 
tistical forecast component of MASS that is unique 
for this application. Section 3 illustrates the perfor- 
mance of MASS for a specific case. Section 4 high- 
lights aspects of a system evaluation including the 
MASS model performance over many cases in pre- 
dicting mesoscale patterns of precipitation. Section 
5 concludes with a summary and a discussion of some 
outstanding issues and problems relating to local me- 
soscale modeling at KSC/CCAS. 

2. System overview 

The MASS forecast system is composed of three 
main components: I ) an initialization module, 2) a dy- 
namical model, and 3) a set of statistical models that 
generate probability forecasts of specific weather 
events from dynamical model output and observa- 
tions. The initialization module and dynamical model 
are summarized in Tables 2 and 3. respectively. 
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Table 2. MASS initialization module attributes. 

Attribute Description Resolution Reference 

Terrain U.S. Central Intelligence Agency (CIA) 5 min — 

global dataset 

Land use/land cover U.l. Geological Survey 30 s Anderson ct aL 

Anderson Level II classification scheme (~1 km) (1976) 


Vegetation Advanced Very High Resolution Radiometer 

data used to compute Normalized Difference 
Vegetation Index 

Analysis based on Antecedent 
Precipitation Index using all available 

days prior to ia4t^^^ion , 



Soil temperature Analysis based on air temperature 

observations typically for 3 or more 
days prior to initialization 2 


1 km Chang and Wetzel 

(1991) 


Variable 
(depending on 
available data) 


Variable 
(depending on 
available data) 


Sea surface temperature Global monthly climatology used as Global dataset 

(SST) first guess for analysis of SST observations (1 degree) 


Objective analysis 


Barnes 

or 

optimum interpolation (OI) 3 


Barnes (1964) 
OI (Gandin 1963) 





virtual temperatures 


1 Presently initialized to a constant value of 0.2. 

: Presently initialized to surface temperature. 

3 Presently used only for sea surface temperature analysis. 


a. Initialization module 

The initialization module, or data preprocessor, 
performs surface parameter specification and surface 
and atmospheric variable initialization. The surface 
parameter routines determine the model horizontal 
grid structure and specify nonprognostic parameters 
such as terrain height, land-water classification, land 
use, and fraction of the surface covered by vegeta- 
tion. The data sources and resolution used to initial- 
ize these parameters are given in Table 2. 

There are a number of in situ and remotely sensed 
data sources that are presently used to initialize the 
MASS model. The gridded data from the National 
Meteorological Center’s (NMC) Nested Grid Model 
(NGM) C grid provides first-guess fields for a Barnes 
(1964) objective analysis of rawinsonde data. The raw 


NGM C-grid data available at KSC/CCAS have a 
horizontal spacing of 1.25° latitude x 2.5° longitude 
on 10 mandatory pressure levels from 1000 to 100 mb. 
MASS incorporates surface data including measure- 
ments of temperature, winds, moisture, and clouds 
from land-based stations, ships, and buoys; and wind, 
temperature, and dewpoint temperature from the 
mesoscale network of instrumented towers surround- 
ing KSC/CCAS. The surface data are objectively 
analyzed to the model grid using a two-pass Barnes 
(1964) objective analysis scheme. The locations 
of available rawinsonde, surface, buoy, ship, and 
KSC/CCAS tower observations at initialization time 
for a typical model run are shown in Fig. 1. 

The three-dimensional initial moisture analyses are 
enhanced by creating synthetic relative humidity 
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Table 3. MASS model attributes. 


Attribute Description Reference 

Boundary-layer physics High-resolution Blackadar Zhang and Anthes (1982) 

Surface energy and moisture Force-restore model; Noilhan and Pianton (1989); 

budget three-layer soil moisture budget equation Mahrt and Pan (1984) 

(cover layer and two soil layers) 

Grid-scale precipitation Diagnostic — condense water vapor — 

physics in excess of supersaturation 1 

Prognostic — conservation equations for cloud Zhang (1989) 

water (ice) and rain water (snow) 
including cloud microphysics 2 

Subgrid-scale precipitation , Kuo type with moist Kuo (1965); Anthes (1977); 

physics convective-scale downdrafts 3 Frank and Cohen (1987) 

; .« ■ 

Fritsch-Chappell with modifications Fritsch and Chappell (1980); 

by Zhang and Fritsch Zhang and Fritsch (1986) 

Radiation physics Free-atmosphere short- and longwave radiation Sasamori (1968); Stephens 

(1978); Savijarvi (1990) 

Initialization None ! 

Lateral boundary conditions Blending with Kreitzberg-Perkey Perkey and Kreitzberg (1976) 

sponge condition 1 
or 

Radiative Orlanski ( 1976) 

Data assimilation Newtonian relaxation Stauffer and Seaman (1990); 

Stauffer etal. (1991) 


1 Used for all 45-km simulations. 

2 Used for all 1 1-km simulations. 

1 Used for all 11- and 45-km simulations. 


(RH) fields from a combination of manually digitized 
radar (MDR) data, visual surface-based cloud obser- 
vations, and infrared satellite data. The scheme con- 
sists of three basic steps and is described in MESO 
(1993) and Young and Zack (1994). 

In the first step, synthetic RH values are derived 
from surface observations of cloud and current 
weather as well as pilot reports of clouds. To obtain 
RH values from surface cloud and weather observa- 
tions, statistical equations that relate visual observa- 
tions of clouds and weather to vertical relative humidity 
profiles were developed from a database of collocated 
surface and rawinsonde observations. An RH-height 
relationship with a vertical resolution of 25 mb was 
derived for each cloud/weather category (e.g., middle 
overcast with precipitation) using the observed cloud- 


base heights as predictors. An objective analysis scheme 
is used to blend these synthetic RH values with RH mea- 
surements using a first-guess gridpoint field of RH. 

The second step uses IR radiance data to estimate 
the fractional cloud coverage and cloud-top height 
distribution in each model grid cell. Cloud base is es- 
timated from the cloud observations at the nearest sur- 
face station. Model grid points are then moistened or 
dried depending on the fractional cloud coverage 
through the use of the same RH-cloud fraction relation- 
ship used to diagnose clouds in the dynamical model. 

In the third step, grid cells with precipitation are 
identified using MDR reports of echo intensity and 
areal coverage of precipitation and the location of 
convective towers determined from the IR satellite 
data following Adler and Negri (1988). The grid cells 
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45 km MASS Grid 



1 1 km MASS Grid 



Fig. 1 . Depiction of the geographical domain covered by the 
horizontal grid matrices used in (a) the 45-km (coarse mesh) and 
(b) the 1 1-km (fine mesh) mesoscale simulations. Panel (b) con- 
tains an expanded view of the 1 1-km domain given by the inner 
rectangle in panel (a). Representative 45- and 1 1-km grid inter- 
vals are labeled in each panel. The locations of available data 
for typical coarse- and fine-grid model runs are shown as solid 
dots for rawinsondes, open squares for surface stations, open dia- 
monds for ships and buoys, and ‘X’s for KSC/CCAS towers. The 


with diagnosed precipitation are brought to near satu- 
ration from the cloud top to the surface of the earth. 

b. Dynamical forecast model 

The dynamical forecast model used in this system 
is version 5.6 of the MASS model. It is a hydrostatic 
three-dimensional primitive equation model that is a 
descendent of version 2.0 described by Kaplan et al. 
(1982). The attributes of the MASS model are sum- 
marized in Table 3. A detailed description of version 
5.6 and specific enhancements to MASS developed 
for application to forecasting at KSC/CCAS are pro- 
vided elsewhere (MESO 1993). MASS 6.0, which 
will include a full nonhydrostatic equation set, a more 
detailed multiphase water microphysics parameteriza- 
tion, and a turbulent kinetic energy-based boundary 
layer parameterization, is scheduled to be completed 
during the second half of 1996. 

The KSC/CCAS real-time version of the model is 
currently run with a coarse-grid spacing of 45 km 
(55 x 50 points) covering the southeastern United 
States and a fine-grid spacing of 1 1 km (45 x 60 
points) covering the Florida peninsula, the eastern 
Gulf of Mexico, and western Atlantic Ocean. There 
may be some concern about executing a hydrostatic 
version of MASS with a grid spacing as fine as 1 1 
km. However, experiments with the nonhydrostatic 
version of the PSU/NCAR model indicate that results 
from nonhydrostatic and hydrostatic simulations are 
virtually identical in most situations for grid spacings 
above 10 km (e.g., Dudhia 1993; Doyle and Warner 
1993b). The extent of the 45- and 1 1-km domains is 
shown in Fig. 1 . The vertical spacing of the model's 
20 sigma layers used for both coarse- and fine-grid 
runs varies from -20 m at the lower boundary (i.e., the 
surface) to - 2 km at the upper boundary (i.e., 100 mb). 

The MASS data preprocessor and model have been 
running twice daily on a four-processor Stardent 3000 
workstation since January 1994. Section 4 provides 
additional details concerning the model’s perfor- 
mance on the Stardent 3000 and other computing plat- 
forms. The attributes and simulation schedule for the 
real-time MASS configuration are summarized in 
Fig. 2. The daily model forecast and data assimila- 
tion schedule consists of two 24-h coarse-grid and 
two 12-h fine-grid runs per day. The 24-h coarse-grid 
run designated COO is initialized with 0000 UTC data 


rawinsonde sites at West Palm Beach (PB1), Tampa Bay (TBW), 
and Cape Canaveral (XMR) used for verification are indicated 
by the three letter station identifiers in panel (b). 
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and assimilates hourly gridded analyses of surface 
and MDR data from 0000 to 0400 UTC. The hourly 
surface analyses used for data assimilation via 
Newtonian relaxation or nudging (Table 3) are de- 
rived from all available synoptic surface, buoy, ship, 
and KSC/CCAS tower observations at the locations 
shown in Fig. 1. The MDR data are transmitted on 
NMC’s Domestic Data Service at 35 minutes past 
each hour. MASS does not presently assimilate any 
asynoptic data available over the coarse- or fine-grid 
domains shown in Fig. 1. The nudging coefficient is 
set to 0.0003 for both surface and MDR analysis 
nudging. Finally, the NGM forecasts generated from 
0000 UTC data are used to derive lateral boundary 
conditions (BC) for the COO run. The BC are linearly 
interpolated in time from the NGM forecast data at 
0, 6, 12, 18, and 24 h. 

The 12-h fine-grid run, designated FI 2, is initial- 
ized with 1200 UTC data and assimilates 1300 UTC 
surface and MDR data. The 1 2-h forecast from COO 
(valid at 1200 UTC) provides the first-guess fields for 
the objective analysis of 1200 UTC data used for FI 2 
initialization. Additionally, the 12-24-h forecast 
fields from COO are used to specify lateral BC for the 
FI 2 run. For each time step of the FI 2 run, the BC 
are linearly interpolated from the COO output at 1-h 
intervals. The cycle is repeated using 1200 UTC 
data to initialize the 24-h coarse-grid run designated 
Cl 2 and 0000 UTC data to initialize the 12-h fine- 
grid run designated F00. 

The main goal of the daily forecast-assimilation 
cycle is to initialize the fine-grid runs as early as pos- 
sible with current upper-air data. Therefore, the F00 
and FI 2 runs are started approximately 1 h after the 
synoptic data times of 0000 and 1 200 UTC, respectively 
(Fig. 2). Since the COO (C 1 2) forecast is designed pri- 
marily to provide first-guess fields and lateral BC for 
the FI 2 (F00) forecast, it is started well after the 
synoptic data time at 0715 UTC (1915 UTC). As a 
result, the 0000 UTC ( 1 200 UTC) NGM initial analy- 
ses and forecasts can be used for the COO (C 1 2) run since 
all of the 0000 UTC (1200 UTC) NGM gridded data 
are usually received by 0300 UTC (1500 UTC) at 
CCAS. The earliest time that forecast products are 
available and the time that all forecast products are 
available from coarse- and fine-grid runs are given 
in Fig. 2. It is important to point out these times are 
for MASS model forecasts executed on an IBM 
RISC 6000/Model 390 rather than the Stardent 3000. 
The same daily forecast and assimilation cycle shown 
in Fig. 2 has been running on the AMU’s Model 390 


Q Daily Forecast / Assimilation Schedule 


Time (UTC) 


~T 

04 


06 


12 

_L 


fcooh 


"I r 

13 16 

Coarse grid (45 km) 


18 

_L 


00 


T 

01 


06 

_L 




Fine grid (11 km) 


EiD" 


Coarse grid (45 km) 




Fine grid (1 1 km) 


■ Data Assimilation Mode 


■Forecast Mode | C00| Cycle Designation 


Approximate Daily Job Schedule ^ 


Time (UTC) 
Cycle start time 


06 

_X_ 


12 

-L 


18 


00 
L 


0105 (F00) 07 1 5 (COO) 1305(F12) 191$(C12) 

I | I | 

Earliest time forecast 0135 (FOO)* 0730(C00)* 1335 (FI2)* 1930(C12)* 
products available ^ I I ^ 

Time all forecast 0300 (F00)* 0830 (C00)* 1500 (F12)* 203ff(C12)* 

products available 


*Note product availability times shown are for MASS model runs on an IBM RISC 6000 / 
Model 390 rather than the Stardent 3000 discussed in the text. 


Fig. 2. Operational real-time daily forecast, data assimilation, 
and job schedule at KSC/CCAS. 


since March 1995. The Model 390 executes MASS 
approximately three times faster than the four-proces- 
sor Stardent 3000. 

c. Statistical model 

The computational constraints and the unavailabil- 
ity of high-resolution initialization data prohibit the 
execution of MASS with sufficient resolution and 
detailed physics to predict precise occurrences of spe- 
cific weather phenomena such as thunderstorms and 
lightning at KSC/CCAS. As a result, a statistical 
model was incorporated into the MASS prediction 
system. The basic concept was to combine model and 
observational data in a way that would permit the 
generation of hourly updates of the probability of 
specific weather phenomena at KSC/CCAS during 
specified time windows. The expectation was that 
model-generated variables would have more predic- 
tive skill in the longer-lead-time forecasts (i.e., early 
in the day) and that the “latest” values of observation- 
based variables would provide most of the informa- 
tion for the short lead-time (a few hours before the 
target time window) forecasts. The system was in- 
tended to provide a mechanism to transition smoothly 
from predictions based more heavily on model-gen- 
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erated variables to those based on observational data 
as the time of the forecast target window approached. 
This approach is similar in concept to the Model 
Output Statistics (MOS) schemes used by NMC to 
generate forecasts of local variables from regional or 
global model output. 

The statistical model consists of a set of linear dis- 
criminant functions (LDFs; Fischer 1938). In the pro- 
totype version of the system, LDFs were developed 
for four consecutive 2-h forecast time windows cov- 
ering the period from 1500 to 2300 UTC and four 
predictand events: 1) a lightning stroke detected 
within 10 km of the KSC/CCAS weather observation 
site (TTS); 2) a report of thunder heard at TTS; 3) a 
report of rain at the TTS site in either regular or spe- 
cial observations; and 4) a report of a wind gust of 
15 m s~' or higher at any of the KSC/CCAS mesonet 
towers within 10 km of TTS. This statistical model 
can be used to generate an estimate of the probabil- 
ity of the occurrence of each event within any of the 
forecast windows. 

The statistical model was designed to use both ob- 
servation-based data and model-generated data simul- 
taneously, generate a new forecast each hour, and 
generate forecasts beginning at 0000 UTC each day 
for the afternoon period (1500-2300 UTC) of that 
day. A separate LDF was constructed for each fore- 
cast-generation hour for each of the predictands. All 
of the selected variables (observation based or model 
generated) that were normally available by the start 
of a particular hour were used as candidate predic- 
tors for that hour. Thus, variables based solely on ob- 
servational data could be included in the prediction 
equation for any hour after the time that they were 
reported. For example, a variable based on the MDR 
data reported at 2035 UTC could be used 25 min af- 
ter the reporting time as a predictor in the 2100 UTC 
forecast equation. In the case of variables computed 
from model-generated data, the variables were eli- 
gible for consideration as a LDF predictor for any 
hour after the time that the model simulation normally 
terminated. Thus, if a scheduled model simulation 
normally began execution at 0230 UTC and finished 
at 0630 UTC, then any variable computed from the 
output of that simulation was considered as a candi- 
date only for the LDFs at or after 0700 UTC. A list 
of the observation-based and model-generated vari- 
ables considered as candidate predictors is given in 
Zack et al. (1993). The predictors for each hour’s LDF 
were selected from the pool of potential predictors by 
evaluating the discriminating power of all combina- 


tions of three variables and selecting the set of three 
that yielded that highest ability to discriminate be- 
tween the occurrence and nonoccurrence of each 
event. The predictor set for each hour was limited to 
three to avoid overfitting of the data in the limited size 
developmental sample. 

A preliminary set of LDFs were derived from a 
sample of 58 warm season cases from the summer of 
1992. The 58 cases were a subset of a sample of 102 
cases for which real-time MASS simulations were 
generated on a daily basis between mid-July and Oc- 
tober of 1992. The sample size for the derivation of 
the statistical equations was set to 58 because that was 
the number of cases for which a complete set of ob- 
servational and simulated data needed to define the 
predictors and predictands was available. The domi- 
nant reason that cases in the 1 02-case database of real- 
time MASS simulations had to be excluded from the 
statistical sample was the inability to retrieve data 
from KSC/CCAS sensors because of communications 
difficulties. As a result, the sample size was undesir- 
ably small. The small sample size prevented the evalu- 
ation of the statistical equations on an independent 
dataset and therefore any results from the statistical com- 
ponent of this system must be viewed as preliminary. 

The observational data and forecast data from the 
45- and 1 1-km simulations were used for the deriva- 
tion of all of the statistical equations. Given the small 
sample size, there was no attempt to isolate the rela- 
tive impact of any data subset (e.g., only 1 1 -km fore- 
casts and observations) on the discriminating power 
of the LDFs. The AMU has compiled simulated and 
observational data from daily real-time MASS runs 
during the warm seasons of 1994 and 1995. This da- 
tabase will permit the statistical equations to be de- 
rived from a larger sample size and will also provide 
an opportunity to evaluate the statistical models on 
an independent data sample. This sample will also 
provide an opportunity to assess the relative value of 
the 45- and 1 1-km simulations in providing predic- 
tors with significant discriminating power for the se- 
lected predictand events. 

An example of the potential impact of the dynami- 
cal-statistical modeling combination on the objective 
forecasting of thunderstorm events at KSC/CCAS is 
illustrated in Fig. 3. This chart illustrates the prob- 
ability of correctly forecasting an event of thunder 
with rain at 1200 UTC in the 2-h period from 2100 
to 2300 UTC during the warm season with four dif- 
ferent methods. The probability estimates are based 
upon the use of decision rules from the LDFs derived 
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from the preliminary sample of 58 warm season cases 
from 1992. The first method is simply an application 
of the climatological probability based solely upon the 
day of the year. The analysis of the 58-case sample 
indicates that this will yield a correct forecast of the 
event about 68% of the time. A more sophisticated 
form of climatological forecast is to combine the cli- 
matology with information from the morning XMR 
sounding. The probability of making a correct fore- 
cast with this method is estimated to be 73%. The third 
method is to use the statistical model with only obser- 
vational data. This approach would be expected to 
generate a correct forecast slightly under 80% of the 
time. The most comprehensive method is to use the 
statistical model with both observation-based and 
mesoscale-model-generated variables. The data from 
the limited 58-case sample indicate that this will give 
the correct forecast slightly over 83% of the time. The 
same performance relationship among these forecasting 
techniques was found to exist for the other time periods 
and events considered in the developmental sample. 


3. Case example 

The 19 February 1992 case provides an illustration 
of the improved forecast guidance that could poten- 
tially be gained by executing a mesoscale model over 
the Florida peninsula. This case was important from an 
operational perspective because the USAF scrubbed 
the second launch attempt of a Delta II rocket from 
Launch Complex 17B at CCAS (see image on cover) 
due to thick clouds (> 4500 ft thick) and disturbed 
weather (i.e., any meteorological phenomena produc- 
ing moderate or greater precipitation). The adverse 
weather was related to an area of thunderstorms that 
developed to the southwest of KSC/CCAS during the 
afternoon hours in advance of a dissipating frontal 
band. The forecasters at CCAS set the overall prob- 
ability of weather constraint violation for the operation 
to 30% just 90 min prior (2029 UTC) to the beginning 
of the launch window. The initial development of this 
isolated convection was not predicted by the NGM 
but was simulated by the MASS model. The perfor- 
mance of MASS for this case was not spectacular, but 
it demonstrates the skill that the model can exhibit 
when mesoscale circulations are an important contribu- 
tor to the initiation and evolution of convective storms. 

At 1200 UTC 19 February 1992, a deep cyclone 
was located over the eastern Great Lakes. A frontal band 
extended southward through coastal South Carolina, 


across the northern portion of the Florida peninsula 
and into the Gulf of Mexico. The band is evident in 
the manually digitized radar (MDR) depictions shown 
in Fig. 4. The band of echoes over northern Florida 
was a result of low to middle clouds and light pre- 
cipitation associated with the frontal zone. This band 
moved very slowly southward and gradually weak- 
ened during the day. The weather to the south of the 
band was generally clear. 

The development of the small area of thunderstorm 
activity to the southwest of KSC/CCAS was appar- 
ently forced by two mesoscale circulations that de- 
veloped over the Florida peninsula during the daylight 
hours. One circulation was associated with a cloud 
boundary on the southern edge of the southwest- 
northeast cloud band over northern Florida. The at- 
mospheric boundary layer was heated significantly in 
the region of nearly clear skies to the south of the 
cloud band while the low-level air within the cloud 
band remained relatively cool. This can be seen in 
Fig. 5 by noting the increase in the surface tempera- 
ture gradient from northern Florida to central Florida 
between 1 500 and 2 1 00 UTC. The observational data 
in Fig. 5 suggest that this north-to-south differential 
heating had a significant impact on the low-level pres- 
sure and wind fields. 

A nested MASS simulation was initialized at 
1200 UTC 19 February 1992. This simulation was ex- 
ecuted over the 45-km coarse-mesh and 1 1 -km fine- 
mesh domains shown in Fig. 1. The model was in the 


100 -A 



Climatology Climatology Statistical Statistical 

(day of year) (with sounding) (observational (observational & 


variables only) model variables) 


Fig. 3. The probability of correctly predicting a thunderstorm 
with rain event for the 2100-2300 UTC at 1 200 UTC using four 
different forecasting techniques. The probabilities are based 
upon the utilization of linear discriminant functions derived from 
a 58-case sample for the summer of 1992. 
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1835 UTC 



Fig. 4. Manually digitized radar depiction for 19 February 1992 at (a) 1535 UTC, 
(b) 1 835 UTC, (c) 2135 UTC, and (d) 2235 UTC. Echo intensities (VIP levels) 1-6 
are shaded gray according to the legend in each panel. 


operational configuration and used only data that are 
routinely available to the operational system. Soil 
moisture can be initialized in MASS from the ob- 
served history of precipitation. However, as in the 
current operational runs at KSC/CCAS, the soil mois- 
ture in this simulation was set to a single climatologi- 
cally representative value throughout the land areas 
in the model domain. Hence, none of the circulation 
features in the simulation can be attributed to varia- 
tions in initial soil moisture. 

Evidence of the circulation associated with the 
cloud boundary can be seen clearly in the model out- 
put data in Fig. 6. The predicted shortwave transmis- 
sivity for 1800 UTC (Fig. 6a) illustrates the simulated 
position of the cloud boundary. At this time, the 
model has already created a substantial north-to-south 
surface thermal gradient (Fig. 6b) in the vicinity of 
the cloud boundary, and there is an incipient meso- 
scale low-pressure center located over central Florida 
on the southern edge of the thermal gradient. The 
model output fields for 2200 UTC (Figs. 6c, d) indi- 
cate that the cloud boundary feature is well developed. 
There is pronounced diffluence in the low-level wind 
field just to the north of the mesolow along the axis 
of a surface pressure ridge in the cloudy region. On 
the northern side of the ridge line, the simulated winds 


are from the southwest, while to the south 
of this ridge the winds have a weaker 
southerly component and a stronger west- 
erly component that results in significant 
confluence farther to the south in the vi- 
cinity of the simulated mesolow. The 
confluence is even stronger in the ob- 
served wind field because the observed 
winds between the ridge and the mesolow 
(Fig. 5d) have a stronger northerly com- 
ponent than the simulated winds (Fig. 6c). 

The second significant mesoscale cir- 
culation was a classic sea breeze. Since 
the large-scale winds on this day were 
from the west, it would be expected that 
the strongest sea-breeze convergence 
would occur along the east coast of the 
peninsula where the sea breeze opposes 
the background large-scale flow. The 
observations at 2200 UTC (Fig. 5d) sug- 
gest that a sea-breeze circulation was 
present along the southeast coast of 
Florida from Miami northward to Vero 
Beach. However, to the north of Vero 
Beach, the sea breeze was suppressed by 
the lack of boundary-layer heating due to cloud cover. 
The simulated 2200 UTC wind field shows an on- 
shore sea-breeze flow present along the southeast 
coast of the peninsula up to the Vero Beach area with 
no onshore flow to the north of the KSC/CCAS area 
as observed (Fig. 5c). 

The small area of new thunderstorm activity that 
was of interest to the KSC/CCAS operational person- 
nel developed between the 2135 and 2235 UTC ra- 
dar summaries (note Figs. 4c, d) near the intersection 
of the sea-breeze convergence line moving westward 
from the east coast and the cloud boundary conver- 
gence line moving southward from the cloud bound- 
ary over north central Florida. The observed position 
of these two convergence zones can be readily in- 
ferred from the surface winds at 2200 UTC (Fig. 5d). 
The model’s moist convective parameterization 
scheme was triggered at approximately 2200 UTC at 
several grid points to the southwest of KSC/CCAS 
and very close to the location where the first radar 
echoes of the new system were observed. No precipi- 
tation was produced by the model anywhere over 
central Florida to the south of the cloud band before 
2200 UTC. The simulated convection moved east- 
ward and crossed the coast just to the south of 
KSC/CCAS by 0000 UTC 20 February as observed. 
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1500 UTC 


1800 UTC 



Fig. 5. Surface weather observation data plotted in standard format at se- 
lected sites and a subjective analysis of altimeter setting (0.5-mb interval) for 
(a) 1500 UTC, (b) 1800 UTC, (c) 2100 UTC, and (d) 2200 UTC. Only a subset 
of the stations used to construct the analysis are shown. 


The simulated 2-h accumulated precipitation for 
the period 2200-0000 UTC (Fig. 6e) indicates the 
area affected by the simulated convective system. The 
maximum simulated precipitation is 1-2 mm just to 
the south of KSC/CCAS, which compares very favor- 
ably with an observation of about 1 mm at Melbourne 
(shown as 0.76 mm in Fig. 7c) and approximately 
2.5 mm at a cooperative observational site just to the 
southwest of Melbourne. The simulated pressure and 
temperature pattern at 0000 UTC (Fig. 6f) indicates 
that the system is accompanied by an area of down- 
draft cooling and a small mesohigh. The timing and 
magnitude of the observed pressure and temperature 
perturbations at Melbourne (not shown) were consis- 
tent with that of the simulated perturbations. 

The 1 2-h precipitation forecasts produced by the 
operational NGM model and the 1 1 -km MASS model 
are compared with observational precipitation data in 
Fig. 7. At the relatively coarse horizontal grid reso- 


lution of 80 km on the C grid, the NGM 
was unable to forecast any of the ob- 
served mesoscale variability of the pre- 
cipitation over the central portion of the 
Florida peninsula. As a result, it drasti- 
cally overpredicted the area covered by 
precipitation. In contrast, the mesoscale 
model was able to forecast a much more 
realistic precipitation distribution. It 
should be noted that current NMC op- 
erational analyses and forecasts (e.g., the 
29-km version of the eta model), which 
were not available when this version of 
the MASS was developed in 1992, em- 
ploy finer horizontal resolution than the 
80-km NGM shown in this comparison. 

This example illustrates a case in 
which the development of moist convec- 
tion was the result of well-defined me- 
soscale features that were attributable to 
differential boundary-layer heating. The 
modeling system tends to perform well 
in this type of scenario since 1 ) many of 
the factors that control the differential 
boundary-layer heating (land-water dis- 
tribution, density of vegetation, soil 
moisture, and cloud patterns) can be rea- 
sonably well mapped for initialization; 
and 2) the heating patterns themselves, 
with the possible exception of those due 
to cloud shading, do not drastically 
change during the course of the simula- 
tion. In contrast, the model does not perform as well 
in cases in which the evolution of convection is 
strongly controlled by the feedback from the convec- 
tion itself (e.g., the development of new convection 
along thunderstorm outflow boundaries). 

4. System evaluation 

While the previous section highlighted the perfor- 
mance of MASS for a specific case, this section fo- 
cuses on the evaluation of the system including 
real-time run statistics, model performance on differ- 
ent computing platforms, and statistical verification 
of the model over many cases. 

cl Real-time MASS run statistics 

The AMU archived output from daily MASS 
model runs from 15 January 1994 through 15 Octo- 
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Fkj. 6. Simulated fields from an 1 1-km MASS simulation initialized at 1200 UTC 19 February 1992: (a) 1800 UTC shortwave 
transmissivity, (b) 1800 UTC 10-m AGL (above ground level) temperatures (dashed lines; 1°C interval) and mean sea level pres- 
sure (solid lines; 0.5 mb interval), (c) 2200 UTC 10-m AGL wind vector and isotachs (1 m s ' interval), (d) same as (b) except for 
2200 UTC, (e) accumulated precipitation (0.5-mm interval) for the 2-h period ending at 0000 UTC 20 February, and (f) same as 
(b) except for 0000 UTC 20 February 1992. 


ber 1994 for the purpose of model evaluation. This 
nine-month archiving period was chosen primarily to 
assess model performance during all four seasons of 
1994. At the end of the archiving period, the number 
of completed MASS model runs was compared with 
the number of total possible runs to measure system 
stability. During this time, no model forecasts were 
lost due to instabilities generated by the model’s phys- 
ics or dynamics, or problems with the model or data 
preprocessor software. Furthermore, the majority of 
45-km runs that were lost were due to hardware prob- 
lems or loss of NGM data used as first-guess fields 
in the MASS preprocessor. In an operational setting. 


MASS would likely be configured to run on a redun- 
dant system and to use alternate first-guess datasets 
such as eta gridded data. In that case, none of these 
45-km forecasts would have been lost. 

Figure 8 displays the percentage of failed runs and 
categorizes the failed runs based on problems with the 
systems (hardware/software) or data. There were a 
total of 462 complete 45-km (coarse grid) runs and 
440 complete 1 1 -km (fine grid) runs out of a total 548 
possible runs. When a coarse-grid run failed, the fine- 
grid run was not executed. At times, the coarse-grid 
run could be restarted and executed at the time that 
the fine grid would normally run. As a result, the 
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Fig. 7. 12-h accumulated precipitation for the period 1200 UTC 19 February-0000 UTC 20 February from (a) the operational 
Nested Grid Model, (b) the MASS 1 1 -km model, and (c) an analysis of observed precipitation from first-order and cooperative 
observing sites. Isopleth interval is 2.54 mm beginning at 0.25 mm. 


number of complete 45- and 1 1-km forecasts do not 
match exactly. 

Figure 8 illustrates that 10.9% of the coarse runs 
were lost due to hardware problems, 2.4% due to soft- 
ware problems, and 2.4% due to loss of data. The 
hardware problems were related to disk and power 
supply failures, while the software problems were 
related to changing the procedures that handle data 
processing. The loss of data includes only NGM 
gridded data that are required as first-guess fields in 
the MASS preprocessor. Figure 8 also shows that of 
the 462 complete 45-km runs, 425 (92%) used NGM 
analysis grids valid at the time of model initialization, 
while 37 (8%) used NGM forecast grids from the pre- 
vious (12-h old) forecast cycle. 

b. Computational performance 

An obvious and critical aspect of any real-time 
numerical forecast system is the requirement for the 
numerical simulations to be generated quickly enough 
to be used as a forecast. Currently, it is not possible 
to simulate satisfactorily, on any computer, all of the 
significant processes that determine the initiation and 
evolution of convection over the entire Florida pen- 
insula in real time. Better resolution and greater de- 
tail in physics are required than can be selected 
for real-time simulations on any workstation or 
supercomputer. The issue is not which computer or 
network of computers will provide the required level 
of computational power, but rather, to what level must 
the simulation model be compromised in order to ex- 


ecute a real-time numerical forecast at a specified 
cost. For example, the PSU, CSU, and FSL real-time 
systems run over very limited domains (Warner and 
Seaman 1990; Cotton etal. 1994; Snook etal. 1995). 
At PSU and CSU, grid nesting is used to concentrate 
the highest horizontal resolution over the local areas 
of interest. In addition, the version of RAMS at CSU 
and FSL is run without using explicit microphysics 
since the scheme greatly increases model run time. 

Greater computational power can always be used 
to improve the quality of the simulations by increas- 
ing the horizontal or vertical resolution and using 
more sophisticated and computationally intensive rep- 
resentations of physical processes. Consequently, the 
capability of a numerical forecast system is strongly 
constrained by the power of the computational plat- 
form that is available for a specific application. How- 
ever, the computational power available for a specific 
price level is rapidly increasing due to continual ad- 
vances in microprocessor and parallel processing 
technology. Hence, a workstation-based numerical 
forecast system should be viewed as a dynamic en- 
tity that improves as processing power improves. 

The cost constraints on the system designed for 
KSC/CCAS dictated that a moderate-cost computer 
workstation be used as the computational platform. 
The selection of the computational platform for the 
current KSC/CCAS system was made in late 1990. 
At that time, the Stardent 3000 workstation was an 
attractive platform for applications, such as MASS, 
that require high-floating point performance. This 
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45 km Runs (548) 



1 1 km Runs (548) 



NGM Gridded Data Usage 
(Complete 45 km Runs) 



37 forecast (8%) 


system operates at a clock speed of 32 MHz and 
is capable of achieving a maximum computational 
rate of 1 floating point operation per clock cycle or 
32 million floating point operations per second per 
CPU board. 

Table 4 illustrates the results of MASS computa- 
tional benchmarks on the 4-CPU Stardent 3000 
system and several current workstation systems as 
well as a Cray 2. Table 4 indicates that the coarse- 
mesh run required just less than 3 h to execute on the 
Stardent 3000 system. However, the time on newer 
computational platforms, which are substantially 
lower in cost, is about one-third of the Stardent 1 s pro- 
cessing time. For example, the 45-km benchmark 
simulation executed in 3284 s (about 55 min) on a 
DEC Alpha Model 3000-600 workstation with a pro- 
cessor clock speed of 175 MHz. This is only 26% 
longer than the execution time on one processor of a 
Cray 2 system. Clearly, the era of low-cost desktop 
supercomputing for numerical weather prediction is 
here. The DEC Alpha Model 3000-600 has a current 
list price of about 25% of the Stardent 3000's 1990 
cost. A similar acceleration in performance is evident 
from the comparison of the 1 1 -km simulations shown 
in Table 4. These results indicate that a similar model 
configuration can be executed much more quickly and 
at a substantially lower cost. Also, a much more sophis- 
ticated model can now be executed in the same time 
and at a lower cost than was possible five years ago. 

c. Statistical verification 

The AMU is completing a detailed evaluation of 
MASS using the nine-month archive (January-Octo- 
ber 1994) of twice-daily runs. The evaluation includes 
objective verification of both gridded and point (or 
station) forecasts, subjective or phenomenological 
verification focusing on selected case studies, and 
rederivation and validation of the LDFs as described 
in section 2c. The results of the complete evaluation will 
be presented in a future paper. This subsection high- 
lights two components of the objective verification. 

The analyses and forecast fields from all available 
coarse-grid, fine-grid, NGM, and persistence fore- 
casts from 15 January 1994 through 15 October 1994 
are bilinearly interpolated to the selected rawinsonde 


Fig. 8. Pie charts showing the percentage of complete 45-km 
runs, complete 1 1-km runs, and the source of first-guess fields 
used for 45-km MASS initialization with analysis denoting the 
0-h NGM gridded analysis data and forecast denoting the 12-h 
NGM gridded forecast data. 
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station locations at West Palm Beach (PBI), Florida; 
Tampa Bay (TBW), Florida; and Cape Canaveral 
(XMR), Florida (see Fig. lb for the location of these 
rawinsonde sites). The average bias and root-mean- 
square errors (rmse) are computed from the twice- 
daily (0000 and 1 200 UTC) rawinsonde observations 
of temperature, relative humidity, wind speed, and 
wind direction at 850, 500, and 300 mb. The errors 
at each pressure level and forecast time are averaged 
for all three stations at both 0000 and 1 200 UTC veri- 
fying times over the entire 9-month period. The rmse 
at 0, 12, and 24 h are given in Table 5. 

The MASS model coarse- and fine-grid analysis 
rmse for temperature and wind speed are typically 
smaller than those from the NGM, indicating that the 
MASS analysis scheme fits the rawinsonde data more 
closely (Table 5). At 12 and 24 h, the errors in the 
NGM and MASS forecasts for temperature, relative 
humidity, and wind speed and direction at 850, 500, 
and 300 mb are similar in magnitude (Table 5). Ad- 
ditionally, an examination of the temperature, wind, 
and moisture bias from the 11- and 45-km MASS 
model forecasts at these same rawinsonde sites (not 
shown) does not reveal any serious systematic errors. 
In general, MASS predicts the large-scale features 
that are sampled by twice-daily rawinsonde observa- 
tions as well as the NGM. Furthermore, the magni- 
tude of the errors for both the NGM and MASS are 
close to the rawinsonde temperature and wind speed 
measurement uncertainty of about 0.6° and 3.1m s" 1 , 
respectively (Schwartz and Benjamin 1995). Thus, it 


would be unrealistic to expect that 
further substantial improvement in 
temperature forecasts could be diag- 
nosed with rawinsonde data. The 
similarity in the error characteristics 
of the two models is not surprising 
since the NGM provides lateral 
boundary conditions for the coarse 
grid, and the coarse grid provides lat- 
eral boundary conditions for fine grid. 
Under strong inflow conditions, the 
information introduced at the lateral 
boundary of the coarse- or fine-grid 
domains can impact the forecasts in 
a relatively short time period. 

The horizontal grid resolution and 
physical parameterizations in MASS 
are likely insufficient to produce 
highly accurate, point-specific fore- 
casts in time or space of warm-sea- 
son convective precipitation. However, to determine 
how well MASS predicts precipitation, both the 
coarse- and fine-grid precipitation forecasts over the 
Florida peninsula are verified using hourly precipi- 
tation data collected by the rain gauge network from 
the St. Johns River, Southwest Florida, and South 
Florida Water Management districts and the gauges 
distributed around KSC/CCAS. The average distance 
between rain gauges is approximately 10 km. The pre- 
cipitation data are analyzed to the 1 1- and 45-km 
model grids using a two-pass Barnes (1964) scheme. 
The hourly gridded precipitation analyses and MASS 
precipitation forecasts are then summed over 1 2 h and 
compared using the equitable threat score (ETS; 
Gandin and Murphy 1992). 

The 0-1 2-h ETS from 1 200 UTC 1 1 -km runs and 
12-24-h ETS from 0000 UTC 45-km runs for each 
precipitation category and month from January 
through October 1994 are shown as bar graphs in 
Fig. 9. [Note that the 0-1 2-h 1 1-km forecasts initial- 
ized at 1200 UTC and 12-24-h 45-km forecasts ini- 
tialized at 0000 UTC are both valid for the same time 
period from 1200 through 0000 UTC (see Fig. 2).] 
The coarse- and fine-grid ETS are less than or equal 
to 0.5 for all thresholds and months (Figs. 9a,b). For 
the most part, fine-grid ETS are larger than coarse- 
grid ETS at thresholds of 2.54, 6.35, and 1 2.7 mm in- 
dicating that the 1 1-km runs are consistently better 
than 45-km runs in predicting precipitation from Janu- 
ary through October 1994. In general, the ETS from 
1 1- and 45-km runs from January through May are 


Table 4. CPU time required for the benchmark simulation of the operational 
configuration of the MASS on different computing platforms. The benchmark 
configuration is essentially identical to the one that is currently used to generate real- 
time simulations at KSC/CCAS. The physics in the coarse- and fine-mesh simulations 
are the same with the exception that the coarse mesh utilizes the diagnostic moisture 
scheme, whereas the fine mesh uses the prognostic moisture scheme. 


45-km coarse mesh: 
24-h simulation 

Computer system (matrix size: 55 x 50 x 20) 

1 1-km fine mesh: 

12-h simulation 
(matrix size: 45 x 60 x 20) 

Cray 2(1 CPU) 

2600 s 

3951 s 

DEC Alpha 3000-600 

3284 s 

5369 s 

IBM RISC 6000/390 

3227 s 

5483 s 

IBM RISC 6000/560 

6293 s 

10 416 s 

Stardent 3000 (4 CPU) 

9718 s 

15 628 s 
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Table 5. Rmse in temperature (°C), relative humidity (%), wind speed (m s and wind direction (degrees) at 300, 500, and 
850 mb for MASS coarse-grid (MASS-C), MASS fine-grid (MASS-F), NGM, and persistence (PERSIS) forecasts. Note that persistence 
errors are computed only at 12 and 24 h, while fine-grid forecast errors are computed only at 0 and 12 h. 


Rmse in temperature (°C) 


Rmse in relative humidity (%) 


Forecast 

hour 

Pressure 
level (mb) 

MASS-C 

MASS-F 

NGM 

PERSIS 

MASS-C 

MASS-F 

NGM 

PERSIS 


300 

0.5 

0.4 

1.2 

— 

12.8 

12.9 

18.4 

— 

0 

500 

0.5 

0.4 

1.1 

— 

12.5 

14.2 

12.8 

— 


850 

0.5 

0.4 

0.9 

— 

9.2 

10.7 

8.4 

— 


300 

1.3 

1.4 

1.2 

1.2 

17.1 

19.8 

26.5 

17.9 

12 

500 

1.2 

1.4 

1.0 

1.1 

19.0 

22.0 

21.1 

22.8 


850 

0.9 

1.0 

: u 

: ’ 1.1 : " 

15.7 

17.8 

14.8 

17.1 


300 

1.5 

— 

1.2 

1.4 

17.7 

— 

30.8 

18.8 

24 

~ 500 

1.4 

— 

Li 

1.4 

20.8 

y. 

22.5 

26.6 


850 

1.0 

— 

1.4 

1.3 

17.7 

— 

16.6 

18.2 


larger than those from June through October (Fig. 9). 
This result suggests that the MASS model provides 
more accurate explicit precipitation forecasts when 
synoptic-scale systems and nonconvective precipita- 
tion dominate the weather in Florida. 

The ETS from 45- and 1 1-km MASS runs are very 
similar to those published for operational models such 
as the NGM and eta model (Junker et al. 1989; 
Zupanski and Mesinger 1995). However, it is impor- 
tant to point out that the skill scores such as the ETS 
do not account for the spatial or temporal errors in 
precipitation forecasts (Olson et al. 1995). For ex- 
ample, the model may predict the correct amount of 
precipitation 2 h later and one grid point farther west 
than observed. In this case, the ETS score would in- 
dicate little or no skill in predicting the event, whereas 
the actual utility of the forecast may be quite good 
considering the spatial and temporal displacement of 
forecast precipitation. 

5. Summary and discussion 

This paper describes the general capabilities and 
operational utility of a version of the MASS that has 
been developed to support operational weather fore- 
casting at KSC/CCAS. For this application, MASS 


has been modified to run on workstation-class com- 
puters. The motivation for running mesoscale mod- 
eling systems like MASS locally at KSC/CCAS is to 
provide detailed short-range (< 24 h) guidance for 
forecasts of winds, clouds, and severe weather such 
as thunderstorms. Forecasting these parameters is im- 
portant for daily operations, and manned and un- 
manned launches and is a challenging task because 
the weather during the warm season at KSC/CCAS 
is dominated by local sea and river breezes, thunder- 
storm outflows, gust fronts, etc. that are forced by 
regional and local factors (i.e., land-water bound- 
aries, vegetation type and amount, soil moisture, etc.). 

The present version of the MASS forecast system 
consists of an initialization module, a dynamical model, 
and a set of statistical models that generate probabil- 
ity forecasts of specific weather events from MASS 
model output and observations. The AMU has been 
running MASS twice-daily on a Stardent 3000 worksta- 
tion since January 1994 and archiving both model out- 
put and observations for the purpose of model evalu- 
ation. During the nine-month period from 15 January 
1994 to 15 October 1994, the largest percentage (10.9%) 
of missed runs resulted from hardware failures. In an 
operational setting, MASS would likely run on a redun- 
dant system, which could have prevented these lost runs. 
Overall, no model forecasts were lost due to instabili- 
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Table 5. ( Continued ) 


Rmse in wind speed (m s _I ) 


Rmse in wind direction (degrees) 


Forecast 

hour 

Pressure 
level (mb) 

MASS-C 

MASS-F 

NGM 

PERSIS 

MASS-C 

MASS-F 

NGM 

PERSIS 


300 

1.4 

1.7 

2.3 

— 

8.5 

13.2 

12.7 

— 

0 

500 

1.1 

1.4 

1.9 

— 

9.8 

17.3 

14.2 

— 


850 

1.0 

1.5 

1.6 

— 

13.0 

26.3 

18.9 

— 


300 

3.3 

3.6 

3.5 

4.5 

19.6 

22.9 

22.4 

28.3 

12 

500 

2.6 

2.9 

2.9 

3.3 

22.9 

23.8 

22.6 

32.8 


850 

2.2 

: : 2J 

2.4 

2.8 

35.2 

41.5 

30.4 

46.4 


300 

3.6 

— 

3.7 

5.8 

22.0 

— 

23.5 

37.2 

; 24 

500 

2.9 

— 

3.2 

4.2 

24.9 

— 

25.1 

36.0 


850 

2.4 

— 

2.6 

3.4 

39.3 

— 

35.8 

53.9 


ties generated by the model’s physics or dynamics or 
problems with the model or data preprocessor soft- 
ware suggesting that MASS is extremely robust and 
would be a very reliable operational system. 

The AMU is completing a detailed evaluation of 
MASS that includes objective verification, case stud- 
ies, and recomputation of the LDFs. A brief exami- 
nation of rmse for temperature, wind, and moisture 
from MASS versus the NGM at selected rawinsonde 
stations over all available cases during the nine-month 
archiving period reveals that MASS is predicting the 
large-scale features as well as the NGM. This result 
is expected since the NGM provides lateral bound- 
ary conditions for the 45-km MASS runs. In fact, veri- 
fication of parameters whose variance is dominated 
by large-scale processes is unlikely to reveal a large 
improvement by mesoscale models such as MASS 
since much of the variance is already accounted for 
by regional-scale models such as the NGM. However, 
to illustrate the utility of running MASS at 1 1 km over 
the Florida peninsula, a case is presented comparing 
the forecast guidance available from MASS versus 
that from the NGM during 19 February 1992 when 
mesoscale circulations were an important contribu- 
tor to the initiation and evolution of convective 
storms. The performance of the 1 1-km MASS run for 
this case, while not spectacular, was superior to the 
80-km NGM forecast especially with respect to the 
distribution of precipitation. 


Although the precipitation forecasts for the one 
case presented here are superior to those from the 
NGM and are reasonable in a qualitative sense, it is 
important to quantify how well MASS predicts pre- 
cipitation from the 45- and 1 1-km runs over many 
cases. For this purpose, the AMU verified precipita- 
tion forecasts from MASS using rain gauge data with 
roughly 10-km spacing over the Florida peninsula. 
The ETS derived from 1200 UTC 1 1-km runs and 
0000 UTC 45-km runs for January through October 
1994 are less than 0.5 and are not consistently better 
than those reported for operational models such as the 
NGM and eta. However, MASS does show greater 
skill as evidenced by higher ETS from January 
through May 1994. It is well known that operational 
models such as the NGM show less skill in forecast- 
ing warm season precipitation associated with small- 
scale convective-type weather systems (Junker et al. 
1989; Olson et al. 1995). 

Based on the ETS derived from mesoscale precipi- 
tation data for MASS model runs, it is apparent that 
precipitation forecasting remains a problem for me- 
soscale models, especially in a subtropical environ- 
ment characterized by weak large-scale forcing such 
as Florida in the warm season. The fact that the 1 1 -km 
MASS runs do not show more skill than operational 
models in forecasting warm season precipitation is 
likely due to a number of factors including insuffi- 
cient horizontal resolution and deficiencies in the 
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MASS runs initialized at 0000 UTC during January through 
October 1994. The ETS for precipitation thresholds of 0.25, 2.54. 
6.35, 12.7. and 25.4 mm from the 45-km coarse-grid forecasts 
are shown in panel (a), while those from the 1 l -km line-grid 
forecasts are shown in panel (b). 


physical parameterizations, especially the Kuo- 
Anthes convective scheme. In addition, the compo- 
nents of the surface energy budget such as evapo- 
transpiration and the representation of the existence 
and impact of subgrid-scale clouds are simplified 
so that MASS can run in real time on workstations. 
Finally, it is difficult to specify accurate mesoscale 
distributions of atmospheric moisture (including 
clouds and preexisting convection), temperature, 
winds, and moisture in the soil and surface cover layer 
from the data sources currently used in MASS. The 
data available from the WSR-88D radars, Doppler 
wind profilers, the new series of geostationary satel- 
lites (GOES-1, J), and Global Positioning System sat- 


ellites (Bevis et al. 1992; Chiswell et al. 1994) may 
offer an opportunity to improve initialization and 
short-range forecasts by MASS if they can be incor- 
porated into the system in real time. 

Workstation-based, real-time numerical modeling 
systems must run fast enough so that the forecasts can 
be used before they expire. However, this obvious and 
critical aspect of these systems must be balanced 
against the desire to improve the quality of the simu- 
lations by increasing the horizontal and/or vertical 
resolution and using more sophisticated physical 
parameterizations. The choice of model attributes 
such as physics, domain size and extent, number of 
grid nests, resolution, and forecast length is highly 
constrained by available computational resources for 
applications of workstation-based real-time modeling 
systems. Clearly, the era of affordable, real-time nu- 
merical weather prediction on workstations has ar- 
rived. Since the monetary cost of computational 
power continues to decrease w ith further advances in 
microprocessor and parallel processing technology, 
there is still opportunity for rapid advancement in 
model performance. Hence, a workstation-based nu- 
merical forecast system should be viewed as a dy- 
namic entity and should evolve in tandem with the 
processing power available at a specilied cost. 

In dosing, if is important to reiterate that MASS 
is not vet an operational system. The AMU has been 
sending real-time MASS output to meteorologists and 
forecasters at the RWO, SMG, and NWS since March 
1 995. As potential end users of model products, their 
feedback to the AMU regarding whether MASS pro- 
vides value-added information to the analysis and 
forecasting of warm season weather at KSC/CCAS 
is a very important component of the overall system 
evaluation. The ultimate recommendation to transi- 
tion MASS for operational use will depend on fore- 
caster feedback and the final results of the AMU\s 
evaluation of MASS. If the recommendation is made 
to transition MASS, there are a number of outstand- 
ing issues involved in the process including documen- 
tation, system certification, training, and lifc-cycle 
costs (i.e., who pays to maintain, repair, and poten- 
tially upgrade the hardware and software). It is an- 
ticipated that the AMU will play an integral role in 
developing a transition strategy that addresses these 
concerns. Furthermore, a transition of MASS into 
operations at KSC/CCAS will require close coordi- 
nation among the AMU, the users, and the organiza- 
tions responsible for certifying and maintaining the 
USAF Eastern Range assets. 
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on the cover 


Delta II rocket with Global Positioning Satellite payload 
on Launch Complex 17B at Cape Canaveral Air Station, 
Florida. The still frame was created from flight line video 
zoomed in on the pad and captures a cloud-to-ground lightning 
strike at 2239:04 UTC on 19 February 1992 while the launch 
count was holding at T-4 minutes. The lightning strike 
occurred offshore about 8 km to the northeast of the pad in 
association with an area of thunderstorms (see case example 
in the article by Manobianco et al. on page 653). 



